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RESEARCH SUMMARY 


A study was conducted on two small watersheds in the Boise 
National Forest to determine the role of forest vegetation in 
maintaining more secure slopes in shallow, coarse-textured 
soils typical of the Idaho batholith. Both soil water piezometry 
and soil shear strength measurements were made in the water- 
sheds. 

Results of the field studies and supporting analyses indicate 
that forest vegetation often provides a critical margin of safety. 
Woody vegetation growing on slopes of the batholith contrib- 
utes to stability by root reinforcement, by soil moisture depletion 
from interception and transpiration, by regulation of snow accu- 
mulation and melt rates, and by soil arching restraint between 
tree stems. Conversely, removal of vegetation from a slope by 
timber harvesting or wildfire results in a loss or reduction of 
effectiveness of these stabilizing mechanisms. Loss of vege- 
tative stabilization in turn can lead to increased frequency of 
landslides as documented in this study. 

Management implications of the study are discussed. Sug- 
gested measures and approaches include more stringent con- 
trols on size and location of clearcut units, greater use of 
“vegetation leave areas” or buffer zones particularly along haul 
roads and next to streams, and construction of hydraulic struc- 
tures that divert water away from critical areas. 
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INTRODUCTION 


Slope Stability Problems 
in the Batholith 


DESCRIPTION OF THE BATHOLITH 

The Idaho batholith is a 16,000 mile? (41 000 km?) moun- 
tainous area of granitic rocks located in central Idaho (see 
location map in fig. 1). Shallow, coarse-textured soils (loamy 
sands to sandy loams) typically develop on slopes that average 
60 percent or more in many drainages. Soils of this type have 
been shown to be extremely erodible (Anderson 1954; André 
and Anderson 1961). 

The granitic bedrock exhibits various degrees of weathering 
and fracturing. An idealized subsurface profile in a slope in the 
natural, undisturbed state, is shown in figure 2. The soil profile 
can be divided roughly into three major zones according to 
degree of weathering, indicated by surface soil (decomposed 
granitics and organic matter), fractured, disintegrated rock, and 
relatively competent, partly weathered bedrock. These zones 
vary in thickness, and changes in composition or texture are 
transitional rather than abrupt. A more precise delineation can 
be made by using a classification for granitic rocks developed 
by Clayton and Arnold (1972). Their classification criteria in- 
clude apparent degree of mineral alteration and rock competen- 
cy with respect to angularity of joint sets and response to a 
hammer blow. 

The surface soil is almost without stones and is derived from 
granitic rock in a fairly advanced state of physical, if not chemi- 
cal, weathering. This horizon usually consists of a sandy, 
coarse-textured soil with little or no intrinsic cohesion. In general 
it would be classified as an A-1-b soil under the American 
Association of State Highway Transportation Officials 
(AASHTO) Classification and as a SW-SM soil under the Uni- 
fied Classification (Gonsior and Gardner 1971). Such relict 
features or vestiges of the original rock structure, as joint planes 
may persist in the soil. It is this part of the profile that is generally 
the most prone or vulnerable to both surficial and mass erosion 
(Durgin 1977). 

Cyclonic storms and/or deep snowpacks that release large 
volumes of water to the soil within short periods are annual 
occurrences in the batholith region. Water supplied by these 
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climatic events rapidly infiltrates the soil surface and continues 
downward until it reaches the zones of reduced hydraulic con- 
ductivity at the interface with the fractured, disintegrated bed- 
rock. Continued inflow of water creates a saturated layer at this 
level, which in turn causes both buildup of a piezometric head 
and subsurface flow along the weathered bedrock surface 
(Megahan 1972). The subsurface flow zone is shown on figure 
2. Depth of the soil piezometric surface is one of the most 
important factors influencing stability of steep mountain slopes 
mantled with shallow, noncohesive soils. 
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Figure 1. — Location map showing the Idaho batholith and study area. 


Figure 2. — Typical soil profile and slope geometry in the Idaho 
batholith (excluding swales, depressions, and other topographic 
lows). Arrows indicate the shallow subsurface flow zone. 


DOCUMENTATION OF SLOPE FAILURES 

The Idaho batholith area contains valuable timber reserves 
as well as other resources important to the regional and nation- 
al economies. This has led to considerable development activ- 
ity, particularly road construction and timber harvesting. The 
batholith is also a critical area with respect to surface and 
subsurface stability. Surface erosion and landslides are preva- 
lent. Consequently, much of this development activity has in- 
creased the occurrence of landslides and provided large ex- 
posed areas from which sediment can readily be removed by 
surface erosion processes (Megahan and Kidd 1972). 

Slope stability problems in the batholith have been studied 
and documented by a number of investigators in the past dec- 
ade (Croft and Adams 1950; Gonsior and Gardner 1971; 
Megahan and Kidd 1972; and Megahan and others 1978). An 
inventory and analysis by Megahan and others (1978) of some 
1,400 landslides on two National Forests reveals some interest- 
ing findings about the nature and cause of these slope failures. 

The survey was carried out on the Clearwater and Boise 
National Forests and represents geologic conditions found in 
the western and central portions of the Northern Rocky Moun- 
tain physiographic province. Most of the area surveyed in the 
landslide inventory falls within the Idaho batholith. Portions of 
the Clearwater National Forest, however, are outside the batho- 
lith. Here, low grade metasediments (mostly quartzite) of the 
Belt Series predominate. 

Almost three quarters (72 percent) of the slides inventoried 
were debris avalanches or debris slides. For a 3-year study 
period, a total of almost 15 acres (6 hectares) of forest land were 
lost to landslides each year on the Clearwater National Forest 
alone. Repair costs to simply clear debris from roads and re- 
place road fill material averaged $56,000 per year. An average 
of 56,000 cubic yards/year (43 000 cubic meters/year) of slide 
material was delivered to active stream channels with adverse 
impacts on the fisheries resources of the region. 

By far the most likely cause of accelerated landslide activity 
was road construction. Roads alone accounted for 58 percent of 
the landslides inventoried. In combination with logging and/or 
forest fires, roads accounted for a total of 88 percent of all 
landslides. Vegetation removal alone accounted for 9 percent of 
the landslide activity. Only 3 percent of the landslides occurred 
on “natural” or undisturbed slopes. These relationships are 
shown graphically in figure 3. 


The Idaho batholith is not unique with regard to its high 
incidence of slope failures. Numerous slides have occurred in 
granitic areas all around the world (Durgin 1977; Jones 1973). 
Granitic masses or batholiths commonly form the core of moun- 
tains that may have extreme topographic relief. Progressive 
physical, chemical, and biological weathering weakens the 
granitic rocks and make them susceptible to erosion and mass 
movement (Clayton and others 1979). Natural events such as 
intense rainstorms, earthquakes, or other perturbations can 
then trigger slides at susceptible sites. Developmental activity 
such as road construction and timber harvesting can accelerate 
or intensify this process. 
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Figure 3.— Landslide occurrence by type of disturbance, western and 
central Northern Rocky Mountain physiographic province (from 
Megahan:and others 1978). 


Potential Effects of Timber Removal 


STABILIZING INFLUENCE 
OF WOODY VEGETATION 

Central to the purpose of this paper is the role of forest 
vegetation in maintaining more secure slopes. This question 
has been analyzed in considerable detail by Gray (1970, 1978), 
who identified four principal mechanisms by which forest 
vegetation enhances stability, namely: 

(a) | Mechanical reinforcement from the rcot system. 

(b) | Regulation of soil moisture content and piezometric 
levels through transpiration, interception, and by 
affecting snow accumulation and rate of melting. 

(c) Buttressing or soil arching action between the trunks 
or stems. 

(d) | Surcharging from the weight of trees. 

The effectiveness of these hydromechanical influences de- 
pends upon soil and slope conditions at particular sites. Root 
reinforcement and buttressing, for example, would be of little 
avail in arresting deep-seated, rotational failures in cohesive 
soils. On the other hand, in shallow, coarse-textured soils which 
are prone to debris sliding and avalanching along an inclined 
bedrock surface, the situation is quite different; root reinforce- 
ment and buttressing may contribute significantly to stability in 
this case. Dramatic examples of such slope stabilization by 
ponderosa pines are shown in figures 4 and 5. The importance 
of these “hydromechanical” contributions of forest vegetation to 
the stability of slopes in the Idaho batholith is examined in 
greater detail later in the paper. 


Figure 4. — Ponderosa pine tree (Pinus ponderosa Laws) buttressing 
and restraining a slope in a road cut in granitic soil of the Idaho 
batholith. Note massive buttress roots in right foreground. 


Figure 5. — Slope buttressing by ponderosa pine. Unbuttressed part of 
slope on left has failed. Mendocino National Forest, California. 


CONSEQUENCES OF REMOVAL 


One way of ascertaining the significance or contribution of 
forest vegetation to the stability of slopes is to examine the 
consequences of its removal. The concensus of investigators 
on this question is clear, namely, indiscriminate removal 
weakens soils and destabilizes slopes. This concensus applies 
not only to the Idaho batholith but to other forested slopes as 
well (Bishop and Stevens 1964; Rice and Krammes 1970; 
Swanston 1974; Swanson and Dyrness 1975; O'Loughlin 1974; 
Swanston and Swanson 1976; and Wu 1976). 


The landslide inventory conducted by Megahan and others 
(1978) in the Boise and Clearwater National Forests revealed 
significant contributions of vegetation roots to stability of steep 
mountain slopes. Landslide hazard was observed to increase in 
direct proportion to the amount of vegetation removed because 
of root decay. The amount of residual vegetation, including both 
trees and shrubs, was an important factor regulating the in- 
crease in landslide hazards following timber cutting or forest 
fire. The rate of root decay relative to the rate of new root growth 
appeared to determine the time of occurrence of maximum 
landslide hazard following vegetation removal. On the average, 
landslide hazards were greatest 4 years after timber removal 
and remained high for about 6 years. By the end of 20 years, 
landslide hazards had returned to their predisturbance levels. 
These slide-vegetation relationships are shown graphically in 
figures 6 and 7. 
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Figure 6. — Landslide occurrence by vegetation crown cover, western 
and central Northern Rocky Mountain province (from Megahan and 
others 1978). 
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Figure 7. — Landslide occurrence by years since logging, western and 
central Northern Rocky Mountain province (from Megahan and 
others 1978). 


DESCRIPTION OF 
PINE CREEK STUDY 


The Pine Creek study was established originally to evaluate 
subsurface flow occurrence on representative granitic slopes 
and the effect of clearcutting and road construction on subsur- 
face flow emerging on the face of road cut slopes. 


Physiographic Setting 


The Pine Creek drainage is a tributary to the Middle Fork of 
the Payette River above Crouch, Idaho. The study was installed 
on two small undisturbed watersheds of 2.4 and 0.8 acres (1.0 
and 0.3 hectares). These microwatersheds are representative 
of first order drainages found in midelevation, fluvial landscapes 
of the Idaho batholitn. 

Slope gradients on the two study watersheds range from 35 
to more than 70 percent and aspects range from northeast to 
northwest (see fig. 8). The soil is a Koppes loamy coarse sand 
and is a member of the sandy-skeletal mixed family of typic 
cryoborolls (Nelson 1976). Sandy loam to loamy sand soils 
overlie moderately weathered and fractured bedrock and range 
in depth from 6 inches (15 cm) on ridges to about 48 inches (122 
cm) in drainage bottoms. Surface soils were almost entirely 
covered by litter up to 1 inch (2.5 centimeters) deep before 
disturbance. Soils are poorly developed, exhibiting only shallow 
A and C horizons. Average physiographic and soils data are 
summarized by watershed in table 1. A general view of the study 
watershed is shown in figure 9. 

Vegetation on the watersheds before disturbance consisted 
of an overstory of ponderosa pine (Pinus ponderosa Laws.), 
Douglas-fir (Pseudotsuga menziesii [Mirb.] Franco), and Engel- 
mann spruce (Picea engelmannii Parry), and an understory of 
small trees and shrubs. 
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Figure 8. — Pine Creek study area showing location of instrument and 
data collection stations in the watersheds. 
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Annual precipitation in the vicinity averages approximately 32 
inches (81 cm). A large portion occurs as snowfalls; conse- 
quently, a large snowpack accumulates often to a total depth in 
excess of 60 inches (152 cm), or 20 inches (51 cm) of water 
equivalent. During the spring, the snowpack melts in a relatively 
short time, contributing large volumes of water to the soil that 
generate subsurface flow and measurable piezometric heads. 


Table 1. — Descriptive data for Pine Creek study watersheds, 
Boise National Forest 


Average 
Average soll 
Watershed slope depth 


Average Average 
Average aspect crown 
elevation (azimuth) cover 


Percent 
Watershed No. 1 
(clearcut) 52 32 


Inches Feet Degrees Percent 


5,050 337NNW_ 63 


Watershed No. 4 
(uncut) 62 30 5,005 335N 43 


Figure 9. — View of Pine Creek study area, Watershed No. 1, Boise 
National Forest, Idaho. Arrows point to slope failures. (Photo taken 
July 1977.) 


History of Land Use 


The study was initiated in 1969 and proceeded on schedule 
until November 1972 when the larger watershed (No. 1) was 
clearcut. A total of 36 ponderosa pines, averaging 25.8 inches 
(65.5 cm) in diameter, and 126 Douglas-fir, averaging 13.4 
inches (34.0 cm) in diameter, were removed as a result of the 
logging. The total volume of timber removal averaged 14,500 
board feet per acre (203 mij*/ha). Slash treatment included 
lopping, scattering, and some hand piling. On August 20, 1973, 
a wildfire swept through the Pine Creek drainage, including the 
two study areas. This unplanned event compromised some of 
the hydrologic studies, but also provided an opportunity to 
compare the postburn response of a cutover and adjacent 
undisturbed watershed. 


Postfire Stability Problems 


Following the wildfire, it was apparent that accelerated sur- 
face soil erosion (by rilling and rain splash) was occurring on the 
clearcut watershed. This was corroborated by studies con- 
ducted by Megahan and Molitor (1975). The erosion was entire- 
ly postfire; there was no evidence of surface erosion on the 
clearcut watershed before the fire. The clearcutting had a two- 
fold effect: it removed protective canopy cover (an important 
source of postfire litter); and it increased surface fuel loading 
that led to a more intense burn on the clearcut watershed. 


Slide 
identification Slide 
number type Length Width 
a aon nnn 2 -n 222-22 ----- Feet------ 

BO-613-75 Debris slide 24 37 
BO-614-75 Debris slide 12 15 
BO-615-75 Debris slide or flow Missing Missing 
BO-616-75 Debris flow(?) 15 15 
BO-617-75 Debris torrent/flow 200 10 
BO-618-75 Debris avalanche 20 10 
BO-619-75 Rotational slump 15 30 
BO-620-75 Rotational slump 35 57 
BO-621-75 Rotational slump and 

debris slide 30 110 
BO-622-75 Rotational slump if 45 
BO-623-75 Debris avalanche 62 10 
BO-624-75 Rotational slump 54 11 
BO-625-75 Debris flow 150 20 
BO-626-75 Debris flow 75 30 
BO-627-75 Debris flow 25 8 
BO-628-75 Debris flow 90 10 
BO-629-75 Debris flow 40 ilZ4 
BO-630-75 Debris flow 65 22 


Mass erosion also occurred in 1974 and 1975 after the fire. It 
was manifested by progressive failures of cut slopes above the 
haul road that traversed the two watersheds (see fig. 9). Other 
landslides also occurred within the 2,200-acre (890-ha) Pine 
Creek burn (table 2). It is interesting that only 50 percent of 
these slides were road associated (most were rotational slumps 
in road cuts). The remainder were shallow, debris slides or flows 
with an average depth of 26 inches (66 cm). These slides 
occurred on natural slopes with an average gradient of 73 
percent (36 degrees). Typical examples of both erosion and 
slope stability problems in the study watersheds are shown in 
figures 10 and 11. 


Table 2. — Landslide study on the Boise National Forest (1975) 


Slides in the Pine Creek Burn 


Depth Volume Gradient Soil depth Cause 
soncnccennennnenen Cubic yard _—_ Percent Inches 
3 50 60 12 Road cut 
es) 7-10 60-70 Missing Not road 
associated 
Missing 9 Missing Missing Road cut 
6 10-12 75 18 Road cut 
1 200 65 16 Not road 
associated 
3 40 Missing Fill failure Road fill 
(culvert) 
1 60 100 22 Road cut 
1 220 80 Missing Road cut 
2 350 70 14 Road cut 
1 30 70 35 Road cut 
1 160 173 Missing Not road 
associated 
2 150 80 Missing Road cut 
2.5 550 80 15 Not road 
associated 
2 240 85 23 Not road 
associated 
3 30 70 16 Not road 
associated 
2 140 70 15 Not road 
associated 
2 55 80 30 Not road 
associated 
2 125 70 16 Not road 


associated 


Figure 10. — Slipout in slope above road cut, Watershed No. 1, Pine 
Creek study area, Boise National Forest, Idaho. (Photo taken July 
1977.) 


Figure 11. — Large gully in slope, Watershed No. 1, Pine Creek study 
area, Boise National Forest, Idaho. Gully was formed by subsur- 
face flow emerging at cut face followed by piping and ‘spring 
sapping,” which triggered localized slope failure. (Photo taken July 
1977.) 


STUDY DESIGN AND METHODS 
Overall Experimental Design 


As noted previously, the original purpose of the Pine Creek 
study was to evaluate subsurface flow occurrence on repre- 
sentative granitic slopes and the effect of clearcutting and road 
construction on subsurface flow. Hydrologic data collection on 
the study watersheds included water inflow, change in water 
storage, and outflow of water occurring as subsurface flow. 

Subsurface flow was measured with collection troughs in- 
stalled on the road at the base of the cut slope across the entire 
width of the watersheds. Estimates of water inflow into the 
watersheds and changes in water storage were obtained from a 
recording rain gage and snow lysimeter. In addition, 45 snow 
stakes and 23 neutron access tubes for soil moisture deter- 
mination were located in a grid pattern on the study watersheds. 
Finally, 25 piezometers were located in suspected water- 
accumulation areas within the watersheds. All data collection 
sites or stations are shown in figure 8. 

Some additional data were collected following the wildfire on 
the study area to evaluate the effects of burning on surface 
erosion rates. Data compilation included rill surveys, a network 
of erosion pins, and collection of eroded material in splash pans 
and in the subsurface flow collection troughs. 

Piezometric levels or ground water conditions are also critical 
determinants of the mass stability of slopes. Accordingly, a 
conjunctive slope stability investigation was undertaken in the 
experimental watersheds in the summer of 1977. Estimates of 
soil shear strength (required for stability analyses) were 
obtained from in-situ borehole shear tests (see fig. 8 for loca- 
tions) and from results of laboratory triaxial tests on granitic soils 
reported by other investigators (Lumb 1962; Gonsior and 
Gardner 1971; Hampton and others 1974": Prellwitz 1975; 
USDA Forest Service 1977). The slope stability investigation 
also included an analysis of the contribution of forest vegetation 
to slope stability from root reinforcement, surcharging, and 
buttressing or soil arching action. 


Piezometer Installation and 
Data Collection 


A total of 25 piezometers were installed in the study water- 
sheds: 14 in Watershed No. 1 and 11 in Watershed No. 4. The 
precise location of these piezometers is shown in figures 12 and 
13. It should be noted that piezometer Nos. 10-14 and 9-11 in 
Watersheds 1 and 4, respectively, were lost in slope failures 
that occurred in 1974 and 1975 (see figs. 12 and 13 for location). 

The piezometers were installed vertically by hand augering a 
2-inch (5-cm) hole through the loamy sand soil to the surface of 
the underlying fractured, disintegrated rock. A typical piezo- 
meter installation is illustrated schematically in figure 14. The 
piezometers were read frequently during the active snowmelt 
period in late spring and at monthly intervals throughout the 
following summer and fall. Both the present and the maximum 
water depths were measured at each piezometer. The max- 
imum reading between measurements was determined by not- 
ing the position of powdered cork on the aluminum rod. The 
powdered cork floats up as piezometric levels rise and adheres 
to the aluminum rod at the point of highest rise. This technique 
was employed by Swanston (1967) in his study of soil water 
piezometry in southeast Alaska. 


‘Hampton, D., W. F. Megahan, and J. L. Clayton. 1974. Soil and rock properties 
research in the Idaho batholith. USDA For. Sci. Lab. Rep., Boise, Idaho 121 p. 
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Figure 13. — Location of piezometers and limits of slide zone, 
Watershed No. 4, Pine Creek study area. 
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Figure 12. — Location of piezometers and limits of slide zone, 
Watershed No. 1, Pine Creek study area. 
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Figure 14. — Schematic illustration of typical piezometer installation, Pine Creek study area. 
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Slope Stability Data Collection 


IN-SITU BOREHOLE SHEAR TESTS 

The borehole, direct shear test is a rapid method for obtaining 
in-situ shear strength parameters. Borehole shear tests avoid 
the need for recovering samples with all the attendant problems 
of sample disturbance and representative sampling (Handy and 
Fox 1967; Wineland 1975). 

The method basically consists of lowering ashear head down 
a 3-inch (8-cm) diameter borehole to a desired depth in the soil 
profile (fig. 15), expanding the shear head against the sides of 
borehole with a Known normal stress, and then recording the 
maximum shear stress required to crank the shear head up the 
hole (fig. 16). The test can be run at several depths in the 
borehole to provide an idea of how the soil shear strength 
parameters vary with depth. 

The boreholes were drilled to the bottom of the soil horizon or 
disintegrated rock interface (fig. 2) using a 2-inch (5-cm) bucket 
auger. The holes were then widened to 3 inches (8 cm) in 
diameter using a thin-walled reamer. Occasionally undisturbed 
samples were taken with the reamer for density and moisture 
content determinations. These samples were also analyzed for 
their grain size distribution. 
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Figure 15. — Expandable shearing head of borehole device prepared 
for lowering down borehole. 
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Figure 16. — Borehole, direct shear test in progress in shallow, coarse- 
textured soil developed on granitic rocks of the Idaho batholith. Pine 
Creek study area, Boise National Forest. 


ROOT DISTRIBUTION AND 
STRENGTH TESTS 

In order to make an assessment of the rooting contribution to 
soil shear strength and slope stability, it was necessary to obtain 
an estimate of both root distribution and root strength in the 
granitic soil of the study watersheds. This estimate was 
obtained from results of studies in granitic soils of other water- 
sheds in the Idaho batholith (Curtis 1964; Burroughs and 
Thomas 1977; and Megahan and others 1978). 

Burroughs and Thomas (1977) measured the concentration 
of intermingled lateral roots (Douglas-fir) in a vertical plane or 
trench midway between trees. Inclusion of all intermingled, 
fresh roots up to 0.4 inch (1 cm) in diameter leads to a calculated 
“root area ratio” of 0.045 percent. Most of the roots counted fell 
into this size class. Roots in the size class 0 to 0.4 inch (0 to 1 
cm) comprised 96 percent of the roots counted in the sample 
area. On the other hand, inclusion of fresh roots up to 3 inches 
(8 cm) in diameter increased the root area ratio to 0.174 percent 
even though these larger roots were not numerous. Root area 
ratios less than 1 percent may not seem significant with regard 
to slope stability, but, as shown later, even a few tenths of a 
percent of root section in the soil is sufficient to provide a 
substantia! and critical amount of shear strength increase in 
shallow, coarse-textured soils. 
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Figure 17. — Tensile resistance of roots of Rocky Mountain Douglas- 
fir as a function of age and size of root (from Burroughs and 
Thomas 1977). 


Burroughs and Thomas (1977) also determined the tensile 
resistance of roots for various size classes and various elapsed 
times after cutting as shown in figure 17. Tensile resistance was 
only measured for Douglas-fir roots in the size interval 0 to 0.4 
inch (0 to 1 cm). The tensile strength computed as tensile 
resistance divided by root cross sectional area (as distinct from 
resistance plotted in fig. 17) decreased from 3,280 Ib/in* (22 600 
kPa) for a0.08-inch (2 mm) diameter root to 2,150 Ib/in* (14 800 
kPa) for a 0.4-inch (10 mm) root. The average tensile strength 
for roots in this size interval, thatis, 0 to 0.4 inch (Oto 1 cm), was 
approximately 2,720 Ib/in* (19 800 kPa). This variation in root 
tensile strength with root size has been observed by other 
investigators as well (Gray 1978; Wu 1976). This type of in- 
formation makes it possible to compute the average tensile 
strength per unit area of soil for a Douglas-fir root-soil system. 
By using asimple root or fiber reinforcement model described in 
the next section, these data can in turn be translated into a 
shear strength increase or pseudo “root cohesion.” 

Unfortunately, Burroughs and Thomas (1977) only looked at 
the concentration of lateral roots. It is the vertical roots (tap and 
sinker roots) that will contribute most to sliding resistance of 
soils on steep, inclined slopes. Studies by Curtis (1964) on 
ponderosa pine indicate deep penetration of tap roots and 
sinker roots in a cylindrical zone around each tree in granitic 
soils of the batholith. The extent of vertical root penetration and 
concentration across the slope surface as a whole is not known 
precisely. 

One way of estimating this concentration is simply to use root 
area ratios that can be calculated from root distributions mea- 
sured by Burroughs and Thomas (1977). Those authors cau- 
tion, however, that their measurements were restricted to lateral 
roots of Douglas-fir crossing a vertical plane midway between 
trees. On the other hand, they also state that the finer roots (0.4 
inch [1 cm] and smaller) are found throughout the rooting zone 
of each tree root system. Furthermore, they note that these 
same roots “.. . will penetrate a shallow soil overlying a weath- 
ered or well-fractured bedrock or glacial till to anchor the soil- 
root mass.” 


More direct data on critical root distribution in granitic soils 
can be obtained from the landslide study by Megahan and 
others (1978). The size and frequency of roots exposed in a 
slide shear plane in a forested slope are summarized in figure 
18. The dominant size class (modal value) of roots in the shear 
plane was 0.20 to 0.40 inches (0.51 to 1.02 cm) in diameter. The 
dominant (modal value) root density class was 6.5 to 10.5 roots 
per square foot (70 to 113 roots per square meter). This in- 
formation yields root area ratios ranging from 0.14 to 0.93 
percent using data from the modal class. Burroughs and Tho- 
mas (1977) found the same class to be dominant for Rocky 
Mountain Douglas-fir in their study, i.e., most of the roots they 
counted fell into the 0 to 0.4 inch (0 to 1 cm) size range. 
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Figure 18. — Landslide occurrence by sheared roots at slide site (from 
Megahan and others 1978). 


In contrast, the root area ratio (0.045 percent) calculated from 
their study for this size class was much lower than those com- 
puted from the data of Megahan and others (1978). The differ- 
ence is most likely caused by the fact that Burroughs and 
Thomas (1977) counted only Douglas-fir roots in vertical planes 
midway between tree stumps, whereas in the other study all 
roots were counted in shear planes parallel to the slope, regard- 
less of species. Wu (1976) measured root area ratios ranging 
from 0.05 to 0.17 and averaging 0.08 percent in his study of 
vegetative influence on landslide occurrence in Southeast Alas- 
ka. His root area ratios were measured at the contact between 
the B and C soil horizons, in soils developed on glacial till slopes 
supporting a spruce-hemlock forest. 

Kozlowski (1971) observed that root structure as well as 
depth and rate of root growth are markedly influenced by the 
rooting medium or environment. The water-holding capacity or 
water availability in a soil is particularly important in this regard. 
Roots tend to avoid regions of high moisture stress and invade 
or permeate moist zones. In the case of shallow, coarse- 
textured soils of the Idaho batholith, the most favorable region 
for roots to exploit from a moisture standpoint is the zone close 
to the contact between the soil and the underlying fractured, 
disintegrated bedrock. In other words, one should expect a fairly 
high concentration of vertical or sinker roots across this contact. 

This expectation is partly confirmed by McMinn (1963) in his 
extensive study of the characteristics of Douglas-fir root sys- 
tems. Hydraulic excavation of root systems in Douglas-fir 
stands revealed a pronounced tendency towards steeply in- 
clined or downward root penetration, not only from the central, 
deep root system, but also from semivertical roots (sinkers) 
from laterals. This trend was most pronounced in older tree 
stands. 


From all the aforementioned studies, a root area ratio in the 
range of 0.05 to 0.15 percent can be considered a reasonable or 
preliminary lower bound estimate in a slope for purposes of 
stability calculations in forested, granitic slopes of the Idaho 
batholith. This topic will be considered further in a subsequent 
section dealing with root reinforcement and its contribution to 
slope stability. 


TREE STEM SURVEY 


Tree stem survey data from nearby experimental watersheds 
in the adjacent Silver Creek drainage were used for calculating 
slope surcharge from the weight of forest vegetation and for 
estimating buttressing or soil restraint from tree trunks firmly 
rooted to the underlying bedrock. 


A stem count and inventory for five cutting units on two 
experimental watersheds in the Silver Creek drainage is 
summarized in table 3. The inventory included only stems 12 
inches (30.5 cm) and over in diameter at breast height. Average 
tree spacings (S) were calculated assuming a simple cubic and 
triangular array. Average size of openings (B) between ‘vertical 
root cylinders” were calculated assuming the diameter of the 
root cyclinder (DR) was approximately five times the stem or 
trunk diameter (D+) at breast height. Field studies by Wu (1976) 
and Curtis (1964) indicate that this diameter ratio is reasonable. 
The various spacing relationships are illustrated schematically 
in figure 19. 


Table 3. — Stem count and tree spacings for experimental watersheds, Silver Creek study area, Boise National Forest 


Number of 
stems 
12 inches 
Total stem in diameter 
Cutting unit Total area area and over 
Acres Feet® 
Control Creek No.3 96.6 2,690 960 
4 5ilkS 3,503 1,567 
5 19.5 1,635 778 
6 2ieil 2,215 880 
K-1 Creek 36.7 1,731 620 


Estimated 
size of 
Estimated? opening 
Average stem Average stem’ root zone between 
diameter spacing diameter root zones 
mons nn enn nena nn nn anc nn nn nne nana nannanasas F@@t ------------------------ 2-2-2 =2 20 nn nono 
1.89 66 9.4 57 
1.69 38 8.4 30 
1.64 33 8.2 25 
1.79 32 8.9 23 
1.88 51 9.4 41 


‘Based on simple, cubic array where spacing (S) = [ (area cutting unit)/(number stems >12 in.) |e 
?Assume Dp = 5D; (estimate based on data in Curtis [1964] and Wu [1976)]). 
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Figure 19. — Schematic diagram of spacing relationships between tree trunks and their ‘vertical root cylinders” in granitic slopes. 
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STABILITY MODELS FOR 
GRANITIC SLOPES 


Infinite Slope Analysis — 
Natural Slopes 


Slopes fail when the shear stress on any potential failure 
surface exceeds the shear strength. It is customary to express 
this balance of forces or stresses in terms of a ratio or factor of 
safety against sliding. This ratio or factor of safety is commonly 
defined as follows: 


F = shear resistance along critical surface 
shear forces promoting sliding on the surface (1) 

Most natural granitic slopes consist of a shallow, cohesion- 
less soil mantle overlying an inclined bedrock contact. An ideal- 
ized soil profile and slope geometry were shown in figure 2. The 
stability of such slopes can be determined by the so-called 
“infinite slope analysis.” This model assumes that the thickness 
of the sliding mass is constant and relatively thin compared to its 
length. This also implies that the sliding surface is parallel to the 
slope and approximately planar over most of its area. The 
infinite slope model has been used previously and with good 
results to analyze the stability of slopes in the Idaho batholith 
(Gonsior and Gardner 1971; Prellwitz 1975). 

According to the infinite slope model, the factor of safety 
against sliding can be expressed by the following mathematical 
equation: 


tand’ 


E +Cpr 
tan B 


tand’'cos’3" (do yi) + (Ya F Hy 


[(Qo + yH) + (sat — Y) Hw] (2) 
F=factor of safety against sliding 
H=thickness (depth) of soil mantle (in vertical direction) 
Hw =height of piezometric surface above bedrock contact 
C,= effective cohesion of soil 
Cr=shear strength increase from root reinforcement 
expressed as a pseudo cohesion 
y=moist density of soil (above piezometric surface) 
Ysat = Saturated density of soil 
Ya = buoyant density of soil (yg = ysat — Yw) 
Yw = density of water 
do=vertical surcharge (from weight of vegetation) 
8 =slope angle or gradient 
o' =effective angle of internal friction of the soil 


Equation 2 is a completely general expression for factor of 
safety, which takes into account the existence or presence of 
cohesion (C,) in the soil, a slope surcharge (q,), and a ground 
water table or piezometric surface in the slope (Hw). In the event 
these are absent, the terms in which they appear are removed 
from the equation. If the soil is completely dry, that is, in the 
absence of any piezometric surface (Hw=O), dry density 
(Yory) May be substituted for moist density (y) in equation 2. 

The relative importance of the various parameters in equation 
2, such as, cohesion as opposed to friction, and the effect on 
safety factor of a change in one variable, such as surcharge, is 
not intuitively obvious. The relative importance of each variable 


inthe equation and direction of change that may be produced by 
altering input variables is best determined by conducting a 
sensitivity analysis using a realistic range of values for each 
input variable. This analysis is conducted for the study water- 
sheds in a subsequent section of the report. 

Slope vegetation and its removal affect several of the input 
variables in equation 2. Most noticeably or obviously affected 
will be the slope surcharge (q,), piezometric height (Hw), and 
root cohesion term (Cr). Reduction in evapo transpiration as a 
result of clearcutting may also affect soil density (y) in addition 
to the piezometric surface elevation. Based on the result of 
soil-root reinforcement studies conducted to date, the angle of 
internal friction of the soil () is affected hardly at all by the 
presence or absence of roots (Gray 1978; Waldron 1977). The 
extent and consequences of forest vegetation and its removal 
on the stability of granitic slopes in the Idaho batholith are 
examined further in the next section of the report. Similar 
studies employing the infinite slope model have been con- 
ducted by Wu (1976), Brown and Sheu (1975), Gray (1978), 
and Ward (1976). 

A simple theoretical model of a fiber-reinforced soil was de- 
veloped by Wu (1976). This model was used by Wu to estimate 
the contribution of rooting strength to slope stability in analyses 
of both forested and cutover slopes in Southeast Alaska. A 
virtually identical model was developed independently by Wal- 
dron (1977) and evaluated in conjunction with direct shear tests 
run in the laboratory on root-permeated homogeneous and 
stratified soil. 

A schematic illustration of the root reinforcement model pro- 
posed by Wu (1976) is shown in figure 20. The model essential- 
ly envisions an elastic root or fiber embedded in a soil matrix and 
initially oriented normal to the shearing surface. Deformation in 
the soil is resisted by tangential forces which develop along the 
fiber and which in turn mobilize the tensile resistance of the 
fiber. These tangential forces (7 on fig. 20) are produced by 
friction or by bonding between the fiber and surrounding soil 
matrix. The soil friction angle () is assumed to be unaffected by 
the reinforcement. This model also assumes that tensile 
strength of the fibers or roots is fully mobilized during failure. 
This requires either fixity of the roots at their ends or roots that 
are long enough and/or frictional enough for the frictional or 


LEGEND 


z = Thickness of shear zone 
x = Horizontal deflection of root 
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Figure 20. — Root reinforcement model. Flexible, elastic root is oriented 
in perpendicular direction to shear surface (after Wu 1976). 


bonding strength between roots and soil matrix to exceed the 
tensile strength of the roots. This assumption appears justified 
from field examination of broken, exposed roots in shear planes 
of landslides (Wu 1976; Megahan and others 1978). 

According to this model the tensile force that develops in the 
fiber can be resolved into a lateral component that directly 
resists shear and into a normal component that increases the 
normal or confining stress on the shear plane. Mathematically 
this translates into an increase in shear resistance or “root 
cohesion” Cr as follows: 


Cr=tr, [cos 8 tan d + sin] (3) 


where Cr=shear strength increase from root 
or fiber reinforcement 
0=angle of shear distortion 
«=angle of internal friction 
ta = average tensile strength of roots per unit 
area of soil 


The average tensile strength of fibers or roots per unit area 
of soil (tz) can be determined by multiplying the tensile strength 
of the roots (Tr) by the fraction of the soil cross section filled or 
occupied by roots (Apr/A). Thus: 


t, =Tr . Ar/A 
Tr=tensile strength of roots, psi 
Apr/A=root area ratio or fraction of soil cross 
sectional area occupied by roots 


(4a) 
where 


This is an extremely useful and important relationship because 
it states that the root contribution to soil strength can be deter- 
mined solely from measuring the tensile strength of the roots, 
Tr, and the fraction of the soil cross section occupied by roots, 
Ap/A. The root cross sectional area, Ar, can be found by count- 
ing the number of roots in different size classes, nj, ina given soil 
cross sectional area, A, and then by summing the product of the 
root numbers in each size class times their corresponding aver- 
age cross sectional area, aj, for that size class. Thus: 
tr =Tr > nj aj (4b) 
A 

where nj=number of roots in size class | 

aj = average cross sectional area of roots in 

size class i, square inches 

A=area of soil in sample count, square inches 

TrR=average tensile strength of roots, psi. 


This linear relationship expressed in equation 4 between root 
tensile strength per unit area of soil and root area ratio has been 
validated by Waldron (1977) for herbaceous plant roots (barley 
and alfalfa) and Wu (1976) for woody plants (spruce and hem- 
lock). 

In the case of natural root systems, the tensile strength tends 
to vary with the size or diameter of the root (Wu 1976; Burroughs 
and Thomas 1977; and Gray 1978). Accordingly, the root ten- 
sile strength term must be included inside the summation and 
the average root tensile strength per unit area of soil com- 
puted by the following relationship: 
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ta= »S T; nj aj (4c) 


A 
where T\=tensile strength of roots in size class i. 

Equations 3 and 4 and the model from which they are derived 
provide a basis for estimating the contribution of roots to soil 
shear strength. The only uncertain or indeterminate variable in 
the equations is the angle of shear distortion (9). This angle will 
vary with the amount of horizontal shear displacement and the 
thickness of the shear zone. 

From results of laboratory direct shear tests conducted by 
Waldron (1977) on various root-permeated soils, the angle of 
shear distortion varied between 40 to 50 degrees. From the 
results of field observations of failures in root-permeated soil 
masses on slopes (Wu 1976), the angle appeared to vary 
between 45 and 70 degrees at most. By running a parametric 
variation or sensitivity analysis on equation 3, Wu (1976) 
showed that the bracketed term is relatively insensitive to all 
expected values of either friction angle () or shear distortion 
angle (9). The bracketed term only varied from 0.92 to 1.31 for 
20 < » < 40 and 40 < © < 70. Thus assuming the midpoint of 
the range to be the most probable value for the bracketed term, 
the shear strength increase from root or fiber reinforcement may 
be estimated to an average or first approximation simply by 


Ca~1.12 tr (Sa) 

or Ca~1.12 TR: Ar (5b) 
A 

or Ca~1.12 = T; nj aj _ (5c) 
A 


This simple theoretical model or relationship permits an esti- 
mate of rooting contribution to shear strength based solely ona 
determination of the root area ratio or concentration of roots ina 
soil cross section and on measurements of tensile strength or 
resistance of the roots themselves. The validity of the model is 
supported by results of direct shear tests run on root-permeated 
soils in both field and laboratory (Endo and Tsuruta 1969; 
Waldron 1977). 

The root reinforcement model can be used to obtain esti- 
mates of root shear-strength increase or ‘root cohesion” (Ca) in 
granitic soils of the Idaho batholith. Root density and tensile 
strength data from the studies by Burroughs and Thomas 
(1977) and Megahan and others (1978) were employed for this 
purpose. Calculated “root cohesions” for Rocky Mountain 
Douglas-fir at various elapsed times after cutting are summa- 
rized in table 4. Root cohesions were calculated using equation 
(5a) and root tensile strength per unit area of soil data for 
Rocky Mountain Douglas-fir reported by Burroughs and 
Thomas (1977). They included data for all roots in the 
size interval 0 to 0.4 inch (0 to 1 cm). As noted previously 
inclusion of roots up to 0.4 inch (1 cm) in diameter corresponds 
to a root area ratio of 0.045 percent. This ratio is a reasonable 
lower bound estimate for root concentration across a potential 
failure surface in the shallow, coarse-textured granitic soils of 
the Idaho batholith. With increasing time after felling, both the 
number of roots and the tensile strength of the remaining roots 
decrease; this accounts for the decrease in root cohesion (table 
4). 
The tabulated data indicate that shear strength increases 
from root reinforcement up to 1.5 Ib/in? (10.3 kPa) are possible in 
granitic soils for initial root area ratios as low as 0.045. Further- 
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more, the data show that some 60 to 70 percent of this strength 
is lost due to root deterioration or decay 5 to 10 years after 
cutting. This elapsed time coincides with the period of greatest 
landslide activity as observed by Megahan and others (1978), 
figure 7. With slightly higher root area ratios (up to 0.15 percent), 
which appear to be entirely possible based on data from Mega- 
han and others (1978) and Wu (1976), the rooting contribution 
to soil shear strength increase will be still higher (assuming 
tensile strength of the roots remains the same). 


Table 4. — ‘‘Root cohesion”’ of soil at various elapsed times after 
felling Douglas-fir for root size class 0 to 1 centimeter 


Residual root cohesion’ (Ca), Ib/inch? 


t = 0 (fresh) t = 1 year t = 5 years t = 10 years 


1.50 0.70 0.50 0.40 


‘Calculated from equation (5a) and root tensile strength data in Burroughs and 
Thomas (1977). 


SOIL ARCHING RESTRAINT MODEL 


Arching in slopes occurs when soil begins to move through 
and around a row of piles (or trees) firmly embedded or 
anchored in an unyielding layer. Under the right conditions, the 
trees are in effect both cantilever piles and abutments of “‘soil 
arches” that form in the ground upslope from the trees. The 
requirement of firm anchoring or embedment of trees in an 
unyielding layer of a slope can occur under the following condi- 
tions: 

(1) Overlying an inclined bedrock contact in shallow residual 
soils or glacial till. 

(2) In sandy slopes, where tree stem bases are deeply 

buried as a result of sand accretion. 
Other conditions pertaining to spacing and diameter of the tree 
trunks, thickness and inclination of the yielding portion of the soil 
profile, and shear strength properties of the soil also determine 
arching effectiveness. 

An arching theory developed for soil slopes by Wang and Yen 
(1974) is based on a semi-infinite slope model and rigid-plastic- 
solid soil behavior. Their theory was developed for a single row 
of embedded piles (of diameter d) spaced a distance (B) apart 
across a slope as shown schematically in figure 21. According 
to this theory, the average arching pressure (p) in a slope and 
the critical spacing (Berit) for arching to occur are given by the 
following equations: 


(mcosBsinB — K,cosBtand — =** cose —mcos’Btand, — os m) 
pee ee ed het ee Se ee ee 


yH 2 K, cos B tan b 


x [ {1 —exp(—2K,ncosBtand)} + 1/2K,exp(—2K,n cosptand) | (6) 


H K, (Ko + 1)tan ob + 2C, 
y 


(7) 


and cet cos (tan B— tan o) — © 


yHcosp 
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Figure 21.— State of plastic deformation and soil arching action around 
a row of piles (trees) embedded in a slope (from Wang and Yen 
1974). 


where B=clear spacing or opening between piles 
Berit = Critical clear spacing between piles 
C,=cohesion in soil 
C;=cohesion along basal sliding surface 
H=depth of yielding or sliding plane 
K,=coefficient of lateral earth pressure at rest 
B/H= relative width or dimensionless spacing 
= X/B=relative distance (in direction of slope) 
p =average lateral pressure or arching pressure 
8=slope angle 
yy = unit weight of soil 
and ,=angle of internal friction in soil and along 
basal sliding surface, respectively. 


= 


The total force (P) developed against a pile of diameter (d) 
embedded in a slope with a thickness or depth of yielding soil 
(H) is given by: 


K K 


Be aod = (yh Bee (8) 


The load on each pile embedded in a slope is thus the summa- 
tion of two loads, one from the pressure at rest of the soil 
immediately uphill from the pile, similar to the lateral pressure 
on a retaining wall. The other is the soil arching pressure trans- 
ferred to the adjacent piles as if each pile is an abutment of an 
arch dam. When the average lateral pressure (p) approaches 
zero, arching action is a maximum. 


RESULTS OF ANALYTICAL AND 
FIELD STUDIES 


Piezometric Responses of Slopes 
to Vegetation Removal 


A number of hydrologic processes are regulated by vegeta- 
tion on a forested slope. For the present study area, important 
processes are snow accumulation and melt and factors affect- 
ing total evaporation losses including interception, transpira- 
tion, and direct evaporation from the soil. These processes in 
turn influence slope stability by helping to regulate the water 
content of the unsaturated soil moisture zone and the depth of 
any underlying piezometric surface. The net effect of vegetation 
removal is to reduce evapotranspiration losses and to increase 
snow accumulation and snowmelt rates thereby creating a 
tendency to raise the piezometric surface. 

Increased soil moisture contents at the end of the growing 
season document the fact that evapotranspiration losses were 
reduced following clearcutting and fire on the study watersheds. 
Also snow survey and snowmelt lysimeter data verify increased 
springtime snow accumulation and snowmelt rates following 
timber removal (Megahan, unpublished data). The combined 
effect of an increased snowpack melting faster in conjunction 
with a greater carryover of fall soil moisture storage causes an 
increase in peak springtime piezometric depths compared to 
the forested condition (fig. 22). The piezometric data shown are 
the annual peak values for the piezometer on each watershed 
that registered the greatest water depth throughout the entire 
study period. The snowpack data are from the nearby Silver 
Creek study area that remained undisturbed throughout the 
course of the present study. 


X Watershed no. 1 (piezometer no. 5) 
@ Watershed no. 4 (piezometer no. 3) 
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Figure 22. — Annual peak piezometric response for peak sample 
points on Watersheds No. 1 and No. 4 and corresponding snow 
accumulations. 


The consistently greater piezometric responses on Water- 
shed No. 1 as compared on Watershed No. 4 reflect the greater 
potential drainage area for piezometers on No. 1 relative to 
No. 4. Water depths on both drainages decreased from 1971 to 
1972 in response to a smaller snowpack in 1972. The 1973 
snowpack was unusually low (there was a 7 percent chance of 
having a snow water content equal to or less than the recorded 
value) and the piezometric levels on the undisturbed Watershed 
No. 4 dropped accordingly. 

A similar trend in piezometer levels did not occur, however, 
on Watershed No. 1, which had been clearcut logged the pre- 
vious November. Rather, the maximum piezometer level 
actually exceeded those recorded in 1972 when snow water 
contents were much greater (there was an 81 percent chance of 
having a snow water content equal to or less than the recorded 
value in 1972). The late fall cutting date for Watershed No. 1 
prevented the development of differences in growing season 
soil moisture storages on the two study watersheds; therefore, 
the increased piezometric levels on Watershed No. 1 are mainly 
the result of differences in snow accumulation and melt during 
the previous winter and spring. 

The following summer, the total soil moisture storage on 
Watershed No. 4 decreased more rapidly than on Watershed 
No. 1 because of reduced evapotranspiration losses caused by 
the timber removal on Watershed No. 1. By mid-August of 1973, 
49 percent more water was stored in the 12- to 60-inch (30- to 
152-cm) depth on Watershed No. 1 (a total of 5.8 inches [14.8 
cm] on Watershed No. 1 and 3.9 inches [9.9 cm] on Watershed 
No. 4). On August 20, 1973, both watersheds were burned in a 
wildfire that killed all nonsprouting vegetation. In the spring of 
1974, piezometric levels were greatly increased on Watershed 
No. 1 relative to predisturbance levels even though the April 1 
snow accumulation was less than the two prelogging years 
studied. The difference reflects both changes in evapotran- 
spiration losses and in snow accumulation and melt rates. 
Piezometric levels also increased on Watershed No. 4, but prob- 
ably not to the maximum because of the development of a 
partial soil moisture deficit prior to the fire. April 1 snow depths 
increased slightly in 1975, relative to 1974, but still were below 
pretreatment amounts. Large increases in piezometric levels 
occurred on both study watersheds in response to combined 
effects of decreased evapotranspiration and snow accumula- 
tion and melt. 

These data do not lend themselves to rigorous statistical 
analysis because of limited sample sizes. A simple graphical 
analysis relating the annual peak piezometer depth to the cor- 
responding April 1 snow water content is informative, however, 
because it suggests the relative effects of vegetation removal 
on piezometric levels (fig. 23). A simple linear regression was 
fitted to the three predisturbance data points for Watershed No. 
4. In spite of the limited data points, the regression had an r 
square value of 0.85 and a standard error of 0.27 feet (0.09 m). 
The regression coefficient was significantly different from zero 
at the 85 percent confidence level. By using the snow water 
contents for 1974 and 1975 to predict the piezometric levels for 
these years and by comparing predicted to measured 
piezometric values, we can estimate increases of about 65 
percent and 100 percent in soil water levels for 1974 and 1975, 
respectively. A similar analysis is not possible for Watershed 
No. 1 because there are only two pretreatment data points. If the 
assumption is made, however, that the regression slope is the 
same as for Watershed No. 4 and that the regression line 
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passes through the average for the two pretreatment years 
(dashed line, fig. 23), we obtain estimates of changes in soil 
moisture levels of 140, 134, and 140 percent for 1973, 1974, 
and 1975, respectively, for Watershed No. 1. 

Estimated increases in maximum piezometric levels range 
from 65 to 140 percent following clearcut timber harvest and/or 
wildfire. We want to emphasize the fact that these increases are 
not statistically significant. The increases, however, are rela- 
tively consistent for both watersheds for all postdisturbance 
years over a wide range of climatic conditions. Moreover, the 
documented onsite hydrologic responses of increased soil 
moisture carryover and increased snow accumulation and melt 
rates following vegetation removal would tend to cause in- 
creased piezometric levels. Everything considered, we feel that 
peak piezometric levels were increased about 100 percent as 
the result of clearcutting and relatively intense wildfire on the 
study area. 


Stability Relationships 


SOIL SHEAR STRENGTH 

The results of borehole shear tests and field density 
measurements in Watershed No. 1 are summarized in table 5. 
Soil friction angles (¢) not only tended to vary with location, but 
with depth as well. Friction angles varied from 29 to 39 degrees, 
and generally increased with depth to a limiting value of 38 to 
39 degrees at the contact between soil (decomposed granitics) 
and fractured, disintegrated bedrock. The friction values re- 
ported in table 5 are station averages for the surface soil alone. 
No cohesion was detected in any of the borehole tests. The 
gradation of a composite sample taken from the 24- to 36-inch 
(61- to 91-cm) depth of the soil horizon in Watershed No. 1 is 
shown in figure 24. The soil consists of 10 percent fine gravel, 
86 percent sand, and 4 percent by weight silt size material: it 
would be classified as a well-graded sand (SW), according to 
the Unified Classification system. In-situ densities ranged from 
92 to 101 Ib/ft* (1.47 to 1.62 g/cm*) with an average value of 96 


ls 


lb/ft? (1.54 g/cm’). Densities tended to increase with depth as did 
the friction angle of the soil. Only the mean density (96 |b/ft®) is 
tabulated in table 5. Subsequent sensitivity analyses showed 
that such small variations in soil density from station to station 
had a negligible influence on calculated slope stability; hence 
the reason for tabulating only the mean. 


Table 5.—Summary of soil-slope data for granitic soil in the Pine 
Creek study watershed 


Station number 


Slope or soil parameters 1 2 3 4 5 6 


Friction angle,’ degrees 34 29 29 34 32 37 
Slope angle, degrees 30 8628 PA eyR Shh 40 
Soil depth, inches 30 «= «36 30) 30, 30) 48 
Soil density,” Ib/ft® 96 96 96 96 96 96 
Dry density,” Ib/ft® 88 88 88 88 88 88 
Saturated density,? Ib/ft® ie ile? SNalre = alaleé alalzé alae 
Cohesion, |b/inch? 0 0 0 0 0 0 


‘In-situ borehole shear test. No cohesion intercept detected. 
?Mean value based on several measurements of field density and water content 
at different stations. 


The soil shear strength parameters, densities, void ratios, 
and gradations measured in the Pine Creek study watersheds 
are comparable to those reported by Gonsior and Gardner 
(1971) during extensive field and laboratory tests of granitic 
soils in the Zena Creek timber sale area in the Payette National 
Forest. Most of the soils investigated by them had slightly higher 
silt size fractions and so were classified as well-graded silty 
sands, or SW-SM materials, according to the Unified Classifica- 
tion. 

Gonsior and Gardner (1971) conducted numerous direct 
shear and triaxial compression tests to investigate the effect of 
unit weight and moisture content upon the shear strength of 
granitic soils. Their test results showed that strengths based 
upon effective stress parameters exhibited negligible variation 
over a wide range of conditions. For design purposes, they 
recommended 35 and 0 as reasonable values for the angle of 
internal friction () and soil cohesion (C,), respectively. Shear 
strength envelopes measured in triaxial compression indicated 
slight residual cohesive strength (on the order of 2 Ib/in® [14 
kPa]) in most tests on saturated specimens. Gonsior and Gard- 
ner (1971) advise, however, that some of the cohesion mea- 
sured in triaxial tests on saturated specimens could be 
accounted for by partial saturation, a curved failure envelope, or 
by membrane resistance; no corrections were made for this 
factor. 

These findings are also corroborated by Lumb (1962) who 
examined the influence of variation in cohesion and friction 
angle with void ratio and degree of saturation. His tests were run 
on both undisturbed and remolded samples of decomposed 
granite with fine, medium, and coarse gradation. Considerable 
“apparent cohesion” was manifest in partially saturated speci- 
mens as a result of capillary forces. As expected, apparent 
cohesion tended towards zero at full saturation (that is, below a 
piezometric surface) in all soils tested, fine or coarse grained. 

Undoubtedly, some residual soil cohesion does exist in grani- 
tic soils; borehole shear test results notwithstanding. In the first 
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Figure 24. — Grain-size distribution of composite soil sample from Watershed No. 1. 


place, the borehole shear test often yields low or nonexistent 
values of cohesion (Wineland 1975). Secondly, full or complete 
saturation is required to eliminate apparent cohesion as a result 
of capillary action. Lastly, some cohesion must be present to 
explain the stability of what otherwise would often be a failed 
slope, particularly in the case of cut slopes. 

Analysis of actual failures in both natural and cut slopes in the 
Idaho batholith by Gonsior and Gardner (1971) suggests that 
cohesion up to 0.9 Ib/in* (6.2 kPa) may be mobilized. A similar 
value was reported by Prellwitz (1975) in his analysis of granitic 
slopes in the batholith. Prellwitz (1975) suggests that a cohe- 
sion of 0.76 Ib/in? (5.2 kPa) is reasonable for soils above the 
phreatic surface and 0.35 Ib/in® (2.4 kPa) below for SW-SM 
materials. 

Much higher cohesions have been reported in laboratory 
triaxial tests On some granitic soils of the batholith (Hampton 
and others 1974).' Sample sites in this study were purposely 
selected to provide a wide range in the weathering properties of 
granitic rocks with the sampling heavier in the more weathered 
rocks. Higher cohesion would be expected under these condi- 
tions because more advanced chemical weathering of the bed- 
rock has occurred, causing formation and accumulation of clay 
colloids. Hydrolysis of mica and feldspars leads to formation of 
clay minerals such as illite and kaolinite. The absence of strong 
solution and eluviation leads to their accumulation. These fac- 
tors may combine to produce a finer grained, more cohesive, 
and correspondingly less frictional soil as shown by the results 
of triaxial compression tests in figure 25. 

Highly weathered bedrock conditions are relatively rare in the 
Idaho batholith because they are generally associated with 
shear zones and zones of secondary hydrothermal alteration. 


Accordingly, high cohesion values such as those shown in 
figure 25 probably occur on less than 5 percent of the upland 
slopes of the batholith. Such soil conditions were not present in 
the Pine Creek study watersheds as evidenced by the shear 
strength test results reported in table 5. 
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Figure 25. — Relationships between angle of internal friction and 
cohesion for various batholith soils (from Hampton and others 
1974, footnote 1). 


SOIL ARCHING RESTRAINT 

Preliminary analyses of field data obtained from forested, 
sandy slopes on the Idaho batholith indicate that these slopes 
meet theoretical criteria for arching restraint between trees. 
Tree spacings, or more importantly the width of openings be- 
tween ‘vertical root cylinders” are of the right order of magni- 
tude for soil arching to manifest itself according to the Wang- 
Yen theory (Wang and Yen 1974). Tree trunks and their associ- 
ated ‘vertical root cylinders,” which are firmly anchored to 
bedrock (fig. 19), potentially can behave as arch abutments. 

Openings between vertical root cylinders appear to average 
around 30 feet (9.1 m) based on stem counts in forested plots on 
the nearby Silver Creek study area (table 3). These data were 
obtained from large survey units which include some unforested 
areas in streams, brush, and rock outcrop. On a smaller more 
localized scale, spacings are considerably less, particularly in 
groves of trees (fig. 26). Based on these last field observations, 
the width of opening between vertical root cylinders in the 
slopes averages about 6 to 7 ft. (1.8 to 2.1 m). 

The maximum allowable opening or critical distance (Bcr) 
between piles (or trees) embedded in a slope can be calculated 
from soil arching theory. This critical distance is shown plotted in 
figure 27 using the soil arching theory for slopes derived by 
Wang and Yen (1974). The critical distance is plotted versus 
cohesion for various assumed values of residual friction and 
cohesion (;, C;) along the basal sliding surface. Other soil and 
slope parameters used in the analysis are typical of shallow 
coarse-textured, granitic soils overlying a steep bedrock rock 
contact (8 = 40°, 6 = 35°,H = 3ft [0.9m] y 100 Ib/ft° 
[1.6 g/cm’}). 


The soil arching analyses show that the critical distance ina 
shallow mantle is very sensitive to cohesion, particularly cohe- 
sion along the basal sliding surface (c,). If no cohesion is 
assumed, and the residual friction (4;) along the basal sliding 
surface is one-half the peak friction (), then the critical spacing 
is 4 ft (1.2 m). On the other hand, if a cohesion (C,) of only 0.35 
lb/in? (2.4 kPa) is assumed with the residual cohesion (c;) along 


Figure 26. — Row of ponderosa pine trees at spacings sufficiently 
close to manifest soil arching restraint between trees. Silver Creek 
study area, Boise National Forest. 
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Figure 27. — Theoretical critical openings (Boca) versus cohesion for 
piles (trees) embedded in a steep, sandy slope (8 = 40°, 6 = 35°, 
H = 3 ft.). Influence of cohesion and friction along the base 
(;, C;) is also Shown. 


the basal sliding surface a mere 12 percent of this value, then 
the critical spacing increases to 21 ft (6.3 m). This distance 
usually approaches or exceeds the size openings between 
“vertical root cylinders” observed in the field, particularly in 
groves of trees (fig. 26). Based on his analysis of granitic slopes 
in the batholith, Prellwitz (1975) suggested that 0.35 Ib/in? (2.4 
kPa) is areasonable lower limit for cohesion of soils beneath the 
phreatic surface. With slightly higher values of cohesion, the 
critical distance increases further thus insuring that soil arching 
effects will be manifest. These values of cohesion in granitic 
soils are well within reason. Possible sources of cohesion or 
apparent cohesion include root reinforcement, cementation, 
clay binder, and capillary stresses (above the phreatic surface). 


SENSITIVITY ANALYSES 

The relative importance of various soil-slope-hydrologic pa- 
rameters on slope stability and the direction of change in slope 
safety factor that may be produced by altering these parameters 
may be determined by conducting sensitivity analyses or para- 
metric variation studies. These same studies also permit eva- 
luation of effect of vegetation removal on slope stability through 
the influence of removal on the parameters themselves. 

There are several types of slope sensitivity analyses that can 
be conducted. All are based in the present case on the infinite 
slope model and on the effect of altering input variables on the 
general stability relationship expressed in equation 2. 

Table 6 includes a summary of stability relationships for each 
of the stations or locations where borehole shear tests were 
conducted. Shown in table 6 are calculated factors of safety for 
each station based on existing or measured values of soil depth 


Table 6. — Summary of stability calculations at different stations 
for granitic soil in the Pine Creek study watershed, 


Idaho 
Station number 
Stability calculations 1 2 3 4 5 6 
Factor of safety’ (F) eee 1045 SOR eA 1.04 0.9 
Critical piezometric level? Glo oat Oy eR ale INA 
(H.,), inch 
Required cohesion for 
stability, Ib/inch* 
H,/H = 0 (dry) 0 0 0) 0 On One 
= 05 0.14 0.24 0.23 0.20 0.20 0.53 
= 1.0 (saturated) 0.40 052 046 0.45 046 0.93 


‘Based on infinite slope model (equation 2.) 
?Piezometric height (above failure surface) at which F = 1.0 for slope and soil 
parameters given in table 5. 


(H), local slope angle (8), soil density (y or yp), friction angle (o), 
and soil cohesion (C,). The calculated factors of safety are close 
to unity and suggest the slope is only marginally secure. The 
calculations are based on zero cohesion (C,; = 0) and dry 
slopes (Hw = 0). Critical piezometric levels (still assuming zero 
cohesion) are also shown. The slope should theoretically fail at 
the critical piezometric level (F = 1.0). These critical 
piezometric levels are on the order of a few inches. The fact that 
the entire slope did not fail when pieometric levels in excess of 
these critical values developed in the slope (figs. 22 and 23) 
means that some cohesion must be present. Required cohe- 
sion to prevent failure at various piezometric levels is also 
calculated and tabulated. Values range from 0.4 to 0.9 Ib/in? (2.8 
to 6.2 kPa) at full saturation (Hy = H). These residual cohe- 
sions are consistent with values reported by Gonsior and Gard- 
ner (1971) and Prellwitz (1975). These observations are based, 
of course, on the assumption that the infinite slope model ade- 
quately represents stability conditions in granitic slopes of the 
Idaho batholith. Limitations of the infinite slope theory in this 
regard are discussed by Hartsog and Martin (1974), but do not 
appear to apply in this case. 

The influence of both friction and cohesion on the factor of 
safety of a typical granitic slope with a shallow soil mantle is 
shown in figures 28 and 29. A slope thickness of 36 inches (92 
cm), slope gradient of 35°, and soil densities of 88 and 117 Ib/ft* 
(1.4 and 1.9 g/cm’), dry and saturated, respectively, were 
selected for the analysis. Factor of safety is plotted against 
cohesion for various values of friction and piezometric elevation 
in the slope. As shown in the figures, stability is far more 
sensitive to soil cohesion than to friction angle, particularly 
when slope becomes fully saturated (Hw = H). It is also clear 
from this analysis that some cohesion must exist in steep slopes 
in order to provide the critical margin for stability when 
piezometric levels rise in the slope. Little cohesion is required to 
maintain a stable slope. Only 0.66 Ib/in? (4.6 kPa) is needed at a 
friction angle of 19° and a slope angle of 35° (tand/tanB = 0.50) 
when Hw/H = 0.5 (fig. 28). For these same conditions of slope 
and friction angle, at full saturation, (Hw/H = 1.0), the required 
cohesion is 0.88 Ib/in? (6.1 kPa) (fig. 29). 

The influence of a rise in piezometric surface on the factor of 
safety (all other factors held constant) can also be determined. 
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Figure 28. — Influence of cohesion on the stability of a sandy residual 
soil resting on an inclined bedrock contact. Hw = 0.5. 
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In this example, a cohesion of 0.5 |b/in? (3.4 kPa) and a friction 
angle of 35° are assumed. When the piezometric surface dou- 
bles in height, the factor of safety drops below one. This 
relationship is shown by the dashed lines in figures 28 and 29, 
respectively. 

Yet another type of sensitivity analysis (Simons and others 
1978)? can be employed that not only reveals the sensitivity, but 
also the direction of change in slope safety factor corresponding 
to a change in any input variable in the infinite slope equation. 
This approach is useful for examining the influence of vege- 
tation removal on slope stability and the influence of changes in 
such variables as surcharge, the effects of which are more 
obscure or counter-intuitive. 

This approach is conducted in four steps. First, a realistic 
range of values (AX;) is selected for each input variable (Xj). 
Second, a base safety factor is computed by using the median 
values for all variables. Third, each input parameter is changed 
across its range of values and a new factor of safety is com- 
puted for each altered input. Fourth, the results are plotted as a 
relative percentage. 

This sensitivity analysis was conducted for the range of input 
variables and their median values shown in table 7. These 
medians are believed typical of conditions for natural slopes in 
forested watersheds in the Idaho batholith. The ranges include 
conditions believed typical of conditions in both forested and 
cutover watersheds. In order to conduct the sensitivity analysis, 
the factor of safety equation was rewritten in a slightly different 
form as follows: 


2(Cs + Cr) leer +QSAT_ 4) oe lore 


F =ywH sin 2B YwH Yw Yw 
Qo is ese \M + Mis (1-M) 
YwH Yw Yw 


where M=relative ground water height (= Fwy 
H 


Table 7. — Range of input variables and their estimated medians 
for soil and slope conditions in watersheds of the 
Idaho batholith 


Input variable Range Median 


Variables not particularly influenced by vegetation: 


B 20-40° 30° 
b 27-42° 35° 
H 12-48 in 30 in 
y 90-120 Ib/ft? 100 Ib/ft® 
YSAT 110-140 lb/ft? 120 lb/ft? 
C, 0-2 Ib/in® 0.75 Ib/in? 


Variables strongly influenced by vegetation: 


do 0-200 Ib/ft? 20 lb/ft? 
Cr 0-1.5 Ib/in® 0.5 Ib/in? 
M 0-1 0.25 


2Simons, D. B., R. M. Li, and T. J. Ward. 1978. Mapping of potential landslide 
areas in terms of slope stability. Report prepared by Eng. Res. Cent., Colo. State 
Univ. for USDA For. Serv., Rocky Mt. For. and Range Exp. Stn., Fort Collins, Colo. 
74 p. 


The results of the sensitivity analysis shown in figure 30 
revealed that some inputs have a linear effect on F while others, 
notably H and B, have strongly nonlinear effects. Factor of 
safety is quite sensitive to both root (Cg) and soil cohesion (Cs). 
In contrast, the slope safety factor is relatively insensitive to 
changes in density (y) and surcharge (loss) q,. Changes in soil 
friction (&) do not have nearly as much influence on safety factor 
as changes in cohesion (Cs and Cr). This finding corrobrates 
the results shown in figures 28 and 29. The influence of relative 
ground water height (M) or piezometric elevation is intermediate 
in effect, except at very high ground water elevations (M——> 1) 
where safety factors decrease sharply. 
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Figure 30. — Percent change in slope safety factor versus percent 
change in input variables. Base safety factor (F) was computed 
using the median or central value for all variables (table 7). 


CONSEQUENCES OF VEGETATION 
REMOVAL 

The preceding sensitivity analyses can be used to help evalu- 
ate the consequences of vegetation removal on the stability of 
slopes in the Idaho batholith. The input variables most strongly 
influenced by vegetation are M, qo, and Cr. Removal of slope 
vegetation tends to decreased root cohesion (Cr), increased 
piezometric levels (M), and decreased slope surcharge (q,). 
The net effect of these changes is to adversely affect stability; 
their extent and significance will be explored further. An excep- 
tion appears to be surcharge that decreases following clearcut- 
ting, which should improve stability based on the sensitivity 
analysis previously discussed. On) the other hand, this im- 
provement is marginal; moreover, it can be shown that at low 
values of cohesion and high ground water elevations surcharge 
has a beneficial influence (Ward 1976). 


impact of Foliage Loss 


Removal of slope vegetation results in a temporary but signifi- 
cant loss of foliage available for interception and transpiration of 
water. This in turn leads to wetter conditions and higher 
piezometric levels in a slope. Results of the Pine Creek study 
support this conclusion (figs. 22 and 23) as do results of other 
hydrologic investigations on effects of timber removal (Gray and 
Brenner 1970; Bethlahmy 1962; Brenner 1973). The impact of 
vegetation removal on soil moisture changes in a slope appears 
to be most critical the first year after cutting. Studies by Hallin 
(1967) showed that, after 3 years, low vegetation that invades a 
cutover site is nearly as effective as old-growth timber in deplet- 
ing moisture. 

The results of soil water piezometry studies in Watersheds 
Nos. 1 and 4 indicate that removal of vegetation by clearcut 
logging can increase piezometric levels as much as 100 per- 
cent. Critical piezometric levels shown in table 6, that is, the 
minimum head required to cause slope failure assuming no soil 
cohesion, were frequently exceeded. The occurrence of slides 
in the slope above the road cut in both watersheds (figs. 10-13) 
and in other watersheds in the vicinity (table 2) reflects the low 
margin of safety under high piezometric conditions. On the 
other hand, the absence of massive and pervasive slope fail- 
ures suggests that some residual cohesion must be present. 
Required cohesions for local stability in Watershed No. 4 for 
different piezometric heads are also shown in table 6. These 
cohesions could be provided in whole or in part by root rein- 
forcement. 


Impact of Root Decay 


The importance of cohesion on stability was clearly estab- 
lished in the preceding sensitivity analyses. The root rein- 
forcement model coupled with root tensile strength and root 
distribution data show that live roots can provide a large fraction 
of the total or apparent cohesion present in granitic soils in the 
batholith. Conversely, studies of root strength loss with time 
after cutting (Burroughs and Thomas 1977) and landslide fre- 
quency with time after cutting (Megahan and others 1978) 
indicate progressive loss of root cohesion following clear- 
cutting. Megahan’s data suggest that landslides are most fre- 
quent 4 to 10 years after logging. These data are consistent with 
findings of other investigators (for example, Bishop and 
Stevens 1964; Swanston and Walkotten 1970). The time of 
minimum stability represents a crossover point between the 
growth and decay curves of root systems of slope vegetation. 
Root strength decline after tree felling is undoubtedly both 
species and site dependent (Burroughs and Thomas 1977). 
The timing or occurrence of slope failures is thus dependent on 
the amount of residual stand on the slope and the rate of 
establishment of new vegetation relative to the root strength 
decline of previously cut trees (Kitamura and Namba 1966). 


Loss of Buttressing and 
Soil Arching Action 


Analysis of spacing relationships and rooting morphology of 
trees in forested slopes of the Idaho batholith indicate the soil 
arching between trees may play an important role in restraining 
soil movement. Several examples of buttressing action by 
embedded tree trunks and root systems were observed (figs. 4 
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and 5). Gonsior and Gardner (1971) reported similar examples 
in their analyses of slope failures in the Idaho batholith. They 
recommended, in fact, that barriers of live trees should remain 
undisturbed immediately below the toe of fill slopes and above 
the cut slope. 

Removal of all large diameter stems by clearcutting, of 
course, gradually eliminates any soil arching restraint or soil 
arching action. The stumps will temporarily provide restraint, 
but when the roots rot and decay these anchor points or “arch 
abutments” will actually become zones of weakness in the 
slope. This will occur because as roots rot and disappear voids 
with no shear strength will be left behind, or infilling with weak 
colloids may occur in the old root channels. In addition, former 
root channels may provide entry points for water and thus 
facilitate rapid buildup of pore pressures. 


Loss of Surcharge 


The sensitivity analysis reported showed that decreasing the 
vertical surcharge (q,.) by removing slope vegetation has a 
beneficial influence on stability, but only a slight one. Under 
certain conditions, surcharge can actually enhance stability. 
Ward (1976) showed that this occurs under the following 
circumstances: 


(C, + Ca)<yw Hw tan  cos’B (10) 


This relationship shows that surcharge is beneficial for low 
cohesion values, high piezometric levels, and relatively gentle 
slopes. Assuming the worst case of maximum rise in 
piezometric surface (Hw = H) and substituting the median 
values in table 6 for the variables on the right-hand side of 
equation (10), a limiting total cohesion of 0.57 Ib/in? (3.9 kPa) 
results. This cohesion is quite possible as an upper, limiting 
value in many granitic slopes. In such slopes, surcharge from 
the weight of trees would have at best a beneficial influence and 
at worst a negligible effect as critical, saturated conditions de- 
velop in the soil. 


MANAGEMENT IMPLICATIONS 


Measures to Minimize 
Mass Erosion Hazard 


LOCATION AND SIZE 
OF CLEARCUT AREAS 

The preceding analyses and findings indicate that many 
slopes in the Idaho batholith are in a state of marginal or 
metastable equilibrium. Such slopes are vulnerable to both 
surficial and mass erosion when vegetation is removed by 
clearcutting or by wildfire. In many instances, road construction 
associated with timber harvesting appears to have a greater 
impact than vegetation removal alone (fig. 3). On the other 
hand, both may have synergistic and cumulative impacts on 
stability that are hard to distinguish and separate. The slope 
failures observed in the slope above the road cut in Watershed 
No. 1 (figs. 10 and 11) are acase in point. The failures appear to 
be associated with the road cut, but may have been caused in 
part by wetter conditions in the slope above and by loss of some 
root cohesion as a result of vegetation removal. 

It is not possible at this point to formulate precise rules for 
location and size of clearcuts to minimize mass erosion 


hazards. The following guidelines are suggested instead: 
1. Limit size of clearcut units; 
2. Stagger location of clearcut units or blocks both in space 
and in time; 
3. Leave buffer zone of trees above and below haul roads; 
4. Leave buffer zone of undisturbed vegetation along all 
streams. 


SELECTION LOGGING 
VERSUS CLEARCUTTING 

The analyses and findings reported here clearly recommend 
leaving as much residual timber stand as possible from the point 
of view of preventing surficial and mass erosion. The greater the 
amount of standing timber, the smaller the amount of soil root 
cohesion loss, the smaller the rise in piezometric levels in a 
slope, and the greater the amount of effective buttressing and 
soil arching action by residual vegetation. All these beneficial 
influences are favored by a selection logging system as 
opposed to clearcutting. 


SITE PREPARATION AND 
ABANDONMENT PROCEDURES 

There are anumber of measures that are routinely employed 
in conjunction with timber harvest operations to minimize slope 
stability problems. These measures are usually specified in 
various State and Federal forest practice rules. They include 
such procedures as the seeding and scarifying of roadbeds, 
removal of temporary road fills, construction of waterbreaks, 
disposal of slash, and establishment of “vegetation leave 
areas.” 

The concept of vegetation leave areas is of particular concern 
and interest in view of the findings reported here. Trees and 
woody vegetation should be left undisturbed in critical areas 
such as steep, slide-prone slopes. Vegetation should also be 
left intact as much as possible along the margins of haul roads 
and streams. Gonsior and Gardner’s (1971) recommendation 
bears repeating in this regard, namely, that barriers of live trees 
should remain undisturbed immediately below the toe of fill 
slopes and above cut slopes. This recommendation should be 
weighed, however, against the likelihood of trees falling across 
roads, owing to possibility of greater vulnerability to root dam- 
age and windthrow after right-of-way-clearing. 

Hydraulic structures should be constructed with regard to 
residual areas of slope vegetation. Crossroad drains and water- 
bars should drain water onto undisturbed vegetation, not over a 
fill slope or into another road or skid trail. Undisturbed vegeta- 
tion should be left to provide water spreading areas large 
enough to accommodate all water draining from roads, skid 
trails, and similar locations. Particular care should be taken to 
avoid “stream piracy” during water spreading operations. This 
can easily happen when water is intercepted by the roadcut in 
one or more microwatersheds and is carried downslope along 
the road in the road drainage system and allowed to spread in 
an adjacent microwatershed. 


General Slope Hazard Rating Scheme 


The recommendation to leave vegetation intact and in-place 
in critical areas during timber harvest operations requires that 
some procedure be employed to identify slopes prone to high 
mass erosion hazard. Several schemes have been devised for 
identifying hazardous slopes (Radbruch-Hall 1976; Ward 1976; 
and Simons and others*). Most of these methods consist of 
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mapping information on slope gradient, soil type, geology, hy- 
drology, and past landslide occurrence. This information is 
integrated by linear combination or factor overlay techniques 
(Hopkins 1977) to produce a composite map of relative slope 
hazard. 

An alternative approach is to base slope hazard ratings on a 
geotechnical model employing principles of limiting equilibrium. 
Geotechnical models such as the infinite slope analysis used 
here (equation 2) explicitly account for the primary factors in 
landslide occurrence such as soil strength, ground water influ- 
ences, vegetative effects, and slope inclination. Geotechnical 
models represent actual field conditions; hence, they can be 
used to analyze the response of a hill slope to temporally and 
spatially varying factors. The geotechnical models or slope 
stability analyses are routinely used by engineers to evaluate 
the stability of a particular hillslope, determine the influence of a 
particular slope modification, and to assess the effectiveness of 
a particular slope protection measure. 

One of the main difficulties with geotechnical models for slope 
hazard analysis is that they are deterministic. As such they do 
not satisfactorily take into account uncertainity and variability in 
the input parameters. A way around this dilemma has been 
developed by Wu (1976), Ward (1976), and Simons and others? 
by casting the stability equation or factor of safety equation ina 
probabilistic framework. Instead of computing a single valued 
safety factor for a slope, one computes a probability of failure. 
Calculated probabilities can then be grouped into three hazard 
classes as suggested by Simons and others’, namely: 

1. High probability when P[F<1] > 60 percent; 

2. Medium probability when 30 <P [F<1] = 60 

percent; 

3. Low probability when P [F<1] < 30 percent. 
where P [F < 1] is the cumulative probability that the safety 
factor (F) is less than or equal to one. 

Computation of the probability of failure requires knowledge 
of the mean and variance of input variables in the safety factor 
equation. This type of information is seldom available without 
extensive testing. This difficulty can be overcome by assuming 
that the input variables are uniformly distributed, random vari- 
ables. With this assumption the mean of a random variable is 
simply found as 


o> ene GO 
2 
and the variance as 


Var [X] = (Xp — Xa)? 
12 

where Xz and Xp are the lower and upper limits on the variable X. 
Thus, probability of failure can be estimated solely from knowl- 
edge of the range in each variable, information which is readily 
available. Simons and others* show that the assumption of a 
uniform distribution provides a conservative estimate of prob- 
ability of failure. The authors also provide a well-documented 
example or application of their method for identifying potential 
landslide areas in terms of their probability of failure. 


CONCLUSIONS 


The following main conclusions can be drawn from the results 
of the Pine Creek study in particular and about the role of forest 
vegetation on stability of slopes in the Idaho batholith in general: 

1. Soils that develop on granitic rocks of the batholith are 
typically shallow, coarse-textured soils (loamy sands to sandy 
loams) that are found on steep slopes that average 60 percent 
or more in many drainages. 

2. Batholith soils tend to be highly erodible and prone to 
mass soil movement particularly when disturbed by road con- 
struction and timber harvesting. 

3. Forest vegetation on the batholith helps to maintain more 
secure slopes by a series of stabilizing mechanisms. These 
include mechanical reinforcement by root systems; soil mois- 
ture depletion by interception, transpiration, and regulation of 
snow accumulation and melting; and by buttressing and soil 
arching action behind embedded tree trunks. 

4. Removal of forest vegetation without regard to slope 
stability can result in loss of the stabilizing influences of forest 
vegetation. Results of the Pine Creek study show that vegeta- 
tion often provides the margin of safety between a secure and 
failed slope. 

5. The factor of safety against sliding in slopes of the batho- 
lith is very sensitive to cohesion. Almost all, or at least a signifi- 
cant fraction of this cohesion, can be provided by root reinforce- 
ment in batholith soils. 

6. Several measures are recommended to mitigate the im- 
pact of vegetation removal on slope stability. These measures 
include selection logging in preference to clearcutting, limitation 
of size of clearcut units, establishment of vegetation leave areas 
in critical areas, and careful integration of diversions and drain- 
age measures with vegetation leave areas. 

7. Live barriers of trees should be left when reasonably 
feasible below the toe of fill slopes and above cut slopes. Buffer 
zones of vegetation should also be ieft along the margin of 
streams. 

8. Critical areas or slopes of high landslide potential can be 
identified by calculating a probability of failure that takes into 
account the uncertainty and variability in the input variables in a 
geotechnical model on which the assessment is based. 
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Results of field studies and supporting theoretical analyses show 
that forest vegetation often provides a significant margin of safety to 
the stability of slopes in the Idaho batholith. Woody vegetation contrib- 
utes to the stability of these slopes by root reinforcement, by soil 
moisture depletion from interception and transpiration, by regulation 
of snow accumulation and melting, and by soil arching restraint be- 
tween tree stems. Removal of vegetation by logging or fire can in- 
crease landslide risks at many locations. 
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The Intermountain Station, headquartered in Ogden, 
Utah, is one of eight regional experiment stations charged 
with providing scientific knowledge to help resource 
Managers meet human needs and protect forest and range 
ecosystems. 

The Intermountain Station includes the States of 
Montana, Idaho, Utah, Nevada, and western Wyoming. 
About 23] million acres, or 85 percent, of the land area in the 
Station territory are classified as forest and rangeland. These 
lands include grasslands, deserts, shrublands, alpine areas, 
and well-stocked forests. They supply fiber for forest in- 
dustries; minerals for energy and industrial development; and 
water for domestic and industrial consumption. They also 
provide recreation opportunities for millions of visitors each 
year. 

Field programs and research work units of the Station 
are maintained in: 


Boise, Idaho 


Bozeman, Montana (in cooperation with Montana 
State University) 


Logan, Utah (in cooperation with Utah State 
University) 


Missoula, Montana (in cooperation with the 
University of Montana) 


a Moscow, Idaho (in cooperation with the Univer- 
a sity of Idaho) 


Provo, Utah (in cooperation with Brigham Young 
University) 


Reno, Nevada (in cooperation with the University 
= Si of Nevada) 


